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Abstract 

We present the results of two-loop calculations of the anomalous dimension matrix for 
the Wilson twist-2 operators in the N=4 Supersymmetric Yang-Mills theory for polarized 
and unpolarized cases. This matrix can be transformed to a triangle form by the same 
similarity transformation as in the leading order. The eigenvalues of the anomalous di- 
^ ■ mension matrix are expressed in terms of an universal function with its argument shifted 
by integer numbers. In the conclusion we discuss relations between the weak and strong 
coupling regimes in the framework of the AdS/CFT correspondence. 
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Parton distributions in QCD satisfy the Balitsky-Fadin-Kuraev-Lipatov (BFKL) pQ 
and Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) [2HH] equations. Next-to-leading 
corrections to the BFKL equation were calculated only recently [3]. It is natural to gen- 
eralize these equations to the supersymmetric case (see Refs. El El and references 
therein). Indeed, the supersymmetric field theories have a number of amazing proper- 
ties, such as a cancellation of quadratic divergencies and non-renormalization theorems 
for interaction terms in the lagrangian. Moreover, the supersymmetry is an excellent 
technical playground for QCD. For example, the empirically established Dokshitzer re- 
lation [2j among elements of the leading order anomalous dimension matrix in the N=l 
supersymmetric limit provides a non-trivial check of results of higher order calculations. 
Another interesting example is the relation between the BFKL and DGLAP equations in 
the N=4 Supersymmetric Yang-Mills (SYM) theory [5JE]- I n this model one can obtain 
the anomalous dimensions of the multiplicatively renormalizable twist-2 operators from 
the eigenvalues of the BFKL kernel (7j. These operators are certain linear combinations 
of the Wilson operators appearing in the theoretical description of the deep-inelastic ep 
scattering El (note, that in the N=4 SYM the beta function is zero and the Bjorken 
scaling for structure functions is strongly violated). Using some assumptions the authors 
of Ref. [0] derived also an expression for the universal anomalous dimension for the N=4 
model in the two-loop approximation. Moreover, the eigenvalues of the anomalous dimen- 
sion matrices in the polarized and unpolarized cases were obtained from this universal 
anomalous dimension by an appropriate integer shift of its argument. In this paper we 
present the results of direct two-loop calculations of these matrices in the N=4 SYM 
theory. 

Now the anomalous dimensions of twist-2 operators in QCD are known up to two loops 
both for the unpolarized [HI El and polarized [TUl E] cases. In the N=4 SYM theory [T2] 
there are one gluon g, four Majorana fermions q, three scalars and three pseudoscalars 
which can be unified in three complex scalars (p. All particles belong to the adjoint 
representation of the gauge group SU(N C ). The transition from QCD to the N=4 SYM 
theory can be performed if one puts in the final expressions Ca = Cf = N c , Tf = 2N C (the 
last substitution follows from the fact, that each gluino qi from four Majorana particles 
gives a half of the Dirac spinor contribution). Furthermore, one should take into account 
the diagrams with virtual scalars in the polarized structure functions and the graphs 
with external scalars in the non-polarized distributions. In the last case the anomalous 
dimension matrix extends to 3 x 3. Below we calculate the anomalous dimensions of the 
following gauge-invariant twist-2 operators: 
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where _D M are covariant derivatives; the spinor \l/ and field tensor G Pfl describe gluinos 
and gluons, respectively, and $ is the complex scalar field appearing in the N=4 super- 
symmetric model. The symbol S implies a symmetrization of the tensor in the Lorenz 
indices fix, and a subtraction of its traces. The anomalous dimension matrices can 
be written as follows for the unpolarized 



7unpol 

and polarized cases 
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Note, that in the super-multiplet of twist- 2 operators there are also operators with fermion 
quantum numbers and operators anti-symmetric in two Lorentz indices . 

Our approach is similar to that of Refs. In particular we calculated unrenormal- 

ized matrix elements of the partonic operators sandwiched between the scalars, fermion 
and gluon states § and the anomalous dimensions were extracted from the expansion of 
the matrix elements through the renormalization group coefficients, with the condition, 
that the renormalized matrix elements satisfy the Callan-Symanzik equations. In our cal- 
culations we used the modified minimal subtraction scheme (MS). Because this scheme 
violates the supersymmetry, the results were transformed to the dimensional reduction 
scheme (DR) JH]) explicitly preserving supersymmetry at least in the two-loop level. For 
this purpose we used the same procedure as in Ref . [Tj|| . Namely, the difference of two- 



loop results in MS and DR-schemes was related to the difference of the finite contributions 
of the corresponding one-loop results. 

In the polarized case one needs an appropriate choice for the 75-prescription. Our 
procedure is analogous to that of Ref. [TU], which based on "reading point" method [20J. 
To begin with, in each trace of the 7-matrix product, we pushed 75 to the right hand side 
using the property of the trace cyclicity. After that we simplified in a straightforward 
way the product of 7-matrices leaving the 75-matrix untouched and used the relation 
1^7^7^7^70-75 = — 4ie pupa . Then the integration over the loop momenta in the space-time 
dimension D = 4 — 2e was performed and the contraction between two Levi-Civita tensors 
(second one appearing from the projector) in four dimensions was done. In the end we 
introduced an additional renormalization constant to restore the anticommutativity of 75 
with other 7-matrices in an accordance with Ref. j21| . 

The final two-loop result for the elements of the anomalous dimension matrix in N=4 
SYM theory in the DR-scheme has the following form (multiplied by a 2 N 2 /(Air) 2 ) in the 

§For the calculations we used the program DIANA |14|. which calls QGRAF ^B] for the generation 
of Feynman diagrams, and the package MINCER [THj for FORM 17 for the evaluation of two-loop 
diagrams. 
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unpolarized case (for even j) 
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the polarized case (for odd j) 
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where 

Q(j) = -I S 1 (j) + 16 S.U) S 2 (j) + 8 S 3 (j) - 8 S 3 (j) + 16 S h2 (j) , (8) 

m = j(^+s2U)+s*u)), (9) 

SuU) = th> ( 10 ) 

i=l 4 

i=l 4 

^0) = Ei^)- (12) 

i=l 

The analytical continuation of functions O'al'' ( J ) ( a ; b — g,q, (p) and 7„^0) (a, b — g,q) to 
the complex values of j can be done analogously to Refs. (221 IE] • The procedure of the 
analytic continuation together with a detailed description of our method of calculations 
will be presented elsewhere. 

The eigenvalues of the anomalous dimension matrices are given below 

7? ) (i) = 7f ) (7) = <Ja-2), (13) 

i^(j) = ii 1] (j) = QU), (14) 

7Skj) = 7 W (j)=Q(j + 2), (15) 
1 §(j) = ^\j) = Q(j-1), (16) 

l { v\j)=l { -\j)=QU + V- (17) 

In fact they coincide with the expressions predicted in Ref. JB]. Indeed, using the two- loop 
result 

7+0') = 7+0' - 1) = 7o(j - 2) = 7-0 - 3) = 7-0 - 4) 



7(j) = -^=ft0-2) + (^)40 



for the universal anomalous dimension 7(j) we can redefine a s — > a s (l — a s N c / (12-k)) to 
remove in Q0) the term proportional to S\{j) (note, however, that it is an additional 
redefinition of a s in comparison with the transition from MS to DR-scheme). After this 
substitution the above universal function Q(j) in two loops coincides with lQQ(j) from 
Ref. 0. 

For the polarized case the Dokshitzer relation is similar to original one (below 7^ 0) = 
7 afe and = % b ) 

%g + If 19 = lm + 2 %i (19) 



and we can find that 

% 9 + \% 9 = QU-i), (20) 

7 S9 -27 OT = QU + 1). (21) 
There are three relations for the unpolarized case 

Igg + Iqg + Isg = Igg + Iqq + Isq = Igs + Iqs + Iss , (22) 

7 OT -47 OT + 3 79S = -^ + 7 9g -^ = ^-^ + 7, s , (23) 

12 lgq - 12 7gs + 3 7(?9 - 3 7gs + 4 7s9 - 4 lsq = (24) 
and one can verify that 

Igg + lig + lsg = Q(j ~ 2) , (25) 

Igg ~ A lgg + 3 7 9 , = Q(j) , (26) 

lgg-lgg- 7 f + 7 f = Q<j + *). (27) 

A complete diagonalization of the above anomalous dimension matrices corresponds to 
the use of a slightly modified basis of the multiplicatively renormalized twist-2 operators 
in comparison with the leading order El due to the breakdown of the superconformal 
invariance (cf. [231 ) - -^ u ^ -^qs. (|19)) - (j27|) are correct in all orders with the replacement of 
Q(j — 2) by the exact universal dimension 

Following the analysis in Ref. jE] the correction to the universal anomalous dimen- 
sions will be constructed as soon as the corresponding QCD anomalous dimensions 
will be calculated (see the recent papers [23| and references therein). 

Recently there was a great progress in the investigation of the N=4 SYM theory 
in a framework of the AdS/CFT correspondence where the strong-coupling limit 
a s N c — > oo is described by a classical supergravity in the anti-de Sitter space AdS^ x S 5 . 
In particular, a very interesting prediction |26| (see also |2Zj) was obtained for the large- j 
behavior of the anomalous dimension for twist-2 operators 

a N 

1 (j) = a(z) In j, z = ^ (28) 

7T 

in the strong coupling regime (see Ref. |2E] for asymptotic corrections): 

lim a = -z 1 ' 2 + ^ + O ( z~ 1/2 ) . (29) 

Here we took into account, that in our normalization 7 (j) contains the extra factor —1/2 
in comparison with that in Ref. |26j . 
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On the other hand, all anomalous dimensions and {i = +, 0, — ) coincide at 
large j and our results for j(J) allow one to find two first terms of the small- ,2 expansion 
of the coefficient a(z) 

71-2 _ I 

lima = -^H z 2 + ... . (30) 

z^o 12 

To go from this expansion to the strong coupling regime we perform a resummation 
of the perturbative result using a method similar to the Pade approximation and taking 
into account, that for large N c the perturbation series has a finite radius of convergency. 
Namely, we present solution of the simple algebraic equation 

- , ~ 1 ~2 /on 

z = -a H — — a . (31) 

From this equation the following large- z behaviour of a is obtained: 

a « -1.1632 z 1/2 + 0.67647 + O (z~ 1/2 \ (32) 

in a rather good agreement with Eq. (j2Uj) based on the AdS / CFT correspondence. Note, 
that if we write for a the more general equation 

2n 

z n = Y J C r a\ (33) 

r=n 

the coefficients C r for n > 3 can be chosen in such way to include all known information 
about a. For n = 2 we can impose on a apart from a correct small- z expansion also the 
condition a = —z 1 ^ 2 for large-z and as a result the asymptotic correction O (1) will be in 
a better agreement with Eq. (|29|) . 

Further, for j — > 2 due to the energy-momentum conservation 

70") = 0'-2)V(2) + ..., (34) 

where the coefficient 7'(2) can be calculated from our results in two first orders of the 
perturbation theory: 

7T 2 

Y(2) = z + 1.2158 z 2 + .... (35) 

6 

Using the same method of resummation as we used above for a, we obtain for large z 

7 '(2) = -0.9071 z 1/2 + 0.6768 + ... . (36) 

Let us take into account, that in this limit 7 = 1/2 + iu + (j — l)/2 — > 1 + (J — 2)/2 
for the principal series of unitary representations of the Mobius group appearing in the 
BFKL equation jl]. Therefore we obtain for large z 

j = 2 - 1.1024 z- l/2 -0.2148 z~ l . (37) 
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in an agreement with the result, that the Pomeron in the strong coupling regime coincides 
with the graviton j^HEO]- The correction ~ z^ 1 ^ 2 to the graviton spin j = 2 coincides in 
form with that obtained in Ref . j2H] from the AdS / CFT correspondence and the coefficient 
in front of z~ x l 2 was not calculated yet. Note, that for the soft Pomeron the correction is 
~ z' 1 [3U]. 

One can attempt to calculate the intercept of the BFKL Pomeron also using its per- 
turbative expansion in Ref. [7j 

3-1 = 2.172Qz - 5.02382 2 . (38) 

After the Pade resummation we obtain in the strong coupling regime j ~ 2.5301 — 
0.84442 -1 in an reasonable agreement with the AdS/CFT estimate (see [30 ). Note, 
however, that in the upper orders of the perturbation theory the BFKL equation should 
be modified by including the contributions from multi-gluon components of the Pomeron 
wave function. 

In the conclusion we want to stress again, that the AdS / CFT correspondence unified 
with a resummation procedure gives a possibility to relate weak and strong coupling 
results. 
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